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The description of the tunnehng of a macroscopic variable in the presence of a bath of locahzed 
spins is a subject of great fundamental and practical interest, and is relevant for many solid-state 
qubit designs. Most of the attention is usually given to the dynamics of the "central spin" (i.e., 
the qubit), while little is known about the spin bath itself. Here we present a detailed study of 
the dynamics of the nuclear spin bath in the Mni2-ac single-molecule magnet, probed by NMR 
experiments down to very low temperatures (T ~ 20 mK) . The results are critically analyzed in the 
framework of the Prokof 'ev-Stamp theory of nuclear-spin mediated quantum tunneling. We find 
that the longitudinal relaxation rate of the ^^Mn nuclei in Mni2-ac becomes roughly T-independent 
below T ~ 0.8 K, and can be strongly suppressed with a longitudinal magnetic field. This is 
consistent with the nuclear relaxation being caused by quantum tunneling of the molecular spin, 
and we attribute the tunneling fluctuations to the minority of fast-relaxing molecules present in the 
sample. The transverse nuclear relaxation is also T-independent for T < 0.8 K, and can be explained 
qualitatively and quantitatively by the dipolar coupling between like nuclei in neighboring molecules. 
This intercluster nuclear spin diffusion mechanism is an essential ingredient for the global relaxation 
of the nuclear spin bath. We also show that the isotopic substitution of by leads to a slower 
nuclear longitudinal relaxation, consistent with the decreased tunneling probability of the molecular 
spin. Finally, we demonstrate that, even at the lowest temperatures - where only T-independent 
quantum tunneling fluctuations are present - the nuclear spins remain in thermal equilibrium with 
the lattice phonons, and we investigate the timescale for their thermal equilibration. After a review 
of the theory of macroscopic spin tunneling in the presence of a spin bath, we argue that most 
of our experimental results are consistent with that theory, but the thermalization of the nuclear 
spins is not. This calls for an extension of the spin bath theory to include the effect of spin-phonon 
couplings in the nuclear-spin mediated tunneling process. 



PACS numbers; 75.45.-|-j, 76.60.-k, 03.65.Yz 
I. INTRODUCTION 



The understanding of quantum tunneling in meso- 
scopic systems has made huge progress in the past 
decades, to the point that nanofabricated devices are now 
being exploited as coherently tunneling two-level systems 
(TLSs) for quantum information purposesJ^i^ Concep- 
tually, a first breakthrough was the proper description 
of the coupling of an effective TLS to an environment 
described by an oscillator bath.^ Whether the system is 
an intrinsic TLS (e.g. a spin s ~ 1/2) or the low-energy 
truncation of a more complicated entity (e.g. the flux 
state of a SQUID), one can generally apply the oscillator 
bath theory when the environment is described by delo- 
calized modes (conduction electrons, phonons, photons, 
etc.) and the couplings of the TLS to each oscillator are 
weak. In many solid-states systems, however, it can be 
necessary to account for localized environmental excita- 



tions whose couplings to the TLS are not weak. This type 
of environment is called "spin bath"— i^!^ and cannot be 
mapped onto an oscillator bath. Importantly, a spin bath 
environment can cause decoherence even at T = and 
is therefore of great relevance for quantum systems that 
are designed to show coherent dynamics, like qubits for 
quantum computation. The prototypical realization of a 
tunneling TLS coupled to a spin bath is the giant spin of 
a single- molecule magnet (SMM)i^i^ii^ These molecular 
systems consist of a core of strongly interacting transition 
metal ions, surrounded by organic ligands. At sufficiently 
low temperatures the core of the molecule behaves effec- 
tively like a single large spin S. When uniaxial magnetic 
anisotropy is present, the reversal of the spin direction 
requires - classically - a large energy, so that the spin di- 
rection can be frozen at very low T. However, in the pres- 
ence of a transverse magnetic field or a biaxial anisotropy, 
the spin direction can be reversed by tunneling through 
the anisotropy barrier pii The electronic spins that form 
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the SMM are magnetically coupled to the nuclear spins 
that either belong to the magnetic ions themselves ('^'^Mn, 
^^Fe, . . . ) or to the surrounding Hgand molecules (^H, 
^^C, . . . ). As a consequence of these couplings, the obser- 
vation of macroscopic quantum tunneling of magnetiza- 
tion in SMMsiSiiSiii^ii^'i^ cannot be understood without 
invoking the dynamics of the nuclear spins themselves.® 
The theoretical predictions for the role of nuclear spins in 
the magnetization tunneling of SMMs^® have been veri- 
fied by a series of experiments on the Fcg compoundj ^^i^^ 
Most remarkably, this material allows to change the iso- 
topic composition of the sample, both by strengthening 
(^^Fe ^^Fe substitution) and weakening (^H 
substitution) the hyperfine couplings, while leaving the 
electronic structure of the SMMs unaffected. As pre- 
dicted, the rate of quantum relaxation of the magneti- 
zation was found to be directly related to the nuclear 
isotopic composition of the sample.— More recently, the 
effect of isotopic substitution has been observed in the 
low-T electronic specific heat of Feg (Ref. [l^ and in 
the dephasing time of coherent electron spin precession 
in CryNi.^" Nuclear spin effects were also invoked in the 
interpretation of /iSR data in isotropic molecules,™ and 
in an alternative description of the short-term magnetic 
relaxation in SMMs.— All these works have analyzed the 
effect of the nuclei on the dynamics of the "central spin" , 
but a crucial aspect of the theory of the spin bath is that 
the tunneling of the central system has repercussion on 
the dynamics of the bath itself, so that the latter cannot 
be simply regarded as an independent source of "noise" . 
Until now, the experiments to probe the electron spin dy- 
namics have not been able to test this delicate aspect of 
the theory. To understand the details of the nuclear spin 
fluctuations, one should then look directly at the nuclear 
spins by means of low-temperature NMR experiments, 
performed under different regimes for the quantum dy- 
namics of the electron spin. These experiments have been 
carried out by several groups , ^'^i^^'^^i^®'^'' but an accurate 
analysis of their implications for the more general theory 
of nuclear-spin mediated quantum tunneling is still lack- 
ing. 

In this work, we present a comprehensive set of ex- 
periments on the dynamics of ^^Mn nuclear spins in the 
Mni2-ac SMM, and we use our results for a critical as- 
sessment of the theory of the spin bath. Our data provide 
definitive proof that the nuclear spin dynamics is strongly 
correlated with that of the central spin, that is, it cannot 
be treated as an independent source of noise. Indeed, we 
find that the nuclear spin fluctuations change dramati- 
cally when the tunneling dynamics of the central spin is 
modified, e.g. by an external magnetic field. In addition, 
we shall demonstrate that the nuclear spins remain in 
thermal equilibrium with the phonon bath down to the 
lowest temperatures (T ~ 20 mK) accessible to our ex- 
periment, where the thermal fluctuations of the electron 
spins are entirely frozen out. This implies that there 
is a mechanism for exchanging energy between nuclei, 
electrons and phonons through the nuclear-spin mediated 



quantum tunneling of the central spin. This is the point 
where the current theoretical description of macroscopic 
quantum tunneling in the presence of a spin bath needs 
to be improved. 

As regards the "macroscopicness" of the quantum ef- 
fects observed in SMMs, we adopt Leggett's view that 
the most stringent criterion is the "disconnectivity" j^^i^^ 
V, which roughly speaking is the number of particles that 
behave differently in the two branches of a quantum su- 
perposition. For instance, while a Cooper pair bo»i is 
a relatively large, lithographically fabricated device, the 
quantum superposition of its charge states involves in fact 
only one Cooper pair, i.e. two electrons, and its discon- 
nectivity is only V = 2. The matter- wave interference in 
fuUerene molecules,'^" for instance, is a much more "quan- 
tum macroscopic" phenomenon, since it means that 60 
X (12 nucleons -I- 6 electrons) — 1080 particles are su- 
perimposed between different paths through a diffraction 
grating. For the spin tunneling in Mni2-ac SMMs dis- 
cussed here, we have 44 electron spins simultaneously 
tunneling between opposite directions, which places this 
system logarithmically halfway between single particles 
and fuUerenes on a macroscopicness scale. 

The paper is organized as follows. Section HIl describes 
the physical properties of the sample used in the exper- 
iments, the design and performance of our measurement 
apparatus, and the methods of data analysis. Section Hill 
presents the experimental results on the nuclear spin dy- 
namics, starting with the NMR spectra, the longitudinal 
and transverse relaxation rates in zero field, and their de- 
pendence on a longitudinal external field. We also study 
the nuclear relaxation in different Mn sites within the 
cluster, and the effect of isotopic substitution in the lig- 
and molecules. In Section IIVI we discuss the thermal 
equilibrium between nuclear spins and phonon bath, the 
experimental challenges in optimizing it, and the indirect 
observation of magnetic avalanches during field sweeps. 
In Section |V] we give an introductory review of the the- 
ory of the spin bath, and apply its predictions to the 
calculation of the nuclear relaxation rate as observed in 
our experiments. Together with the information on the 
thermal equilibrium of the nuclear spins, this will allow 
us to draw clear-cut conclusions on the status of our cur- 
rent theoretical understanding of quantum tunneling of 
magnetization. We conclude with a summary and impli- 
cations of the results in Section IVll 



II. EXPERIMENT 

A. Sample properties 

We chose to focus our study on the well- 
known [Mni20i2(02CMe)i6(H20)4] (Mni2-ac) com- 
pound, which belongs to the family of SMMs with the 
highest anisotropy barrier. As we shall see below, the 
rationale for choosing a SMM with high anisotropy bar- 
rier is that the electron spin fluctuations become slow 
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FIG. 1: (Color online) (a) Structure of the Mni2-ac duster, with the labelling of the three inequivalent Mn sites as described in 
the text, (b) Energy level scheme for the electron spin as obtained from the Hamiltonian ([T}, retaining only the terms diagonal 
in Sz ■ The non-diagonal terms allow transitions between states on opposite sides of the anisotropy barrier by means of quantum 
tunneling (QT). In the presence of intrawell transitions induced by spin-phonon interaction (S-Ph), thermally assisted quantum 
tunneling (Th-A T) between excited doublets can also take place. 



on the NMR timescale already at temperatures of a few 
kelvin. The structure of the cluster'^^ (Fig. [l]) consists of 
a core of 4 Mn'*+ ions with electron spin s = 3/2, which 
we shall denote as Mn^^^, and 8 Mn^+ ions (s = 2) on 
two inequivalent crystallographic sites, Mn'^' and Mn'^' 
[Fig. [DJa)]. Within the molecular cluster, the electron 
spins are coupled by mutual superexchange interactions, 
the strongest being the antiferromagnetic interaction be- 
tween Mn(i) and Mn^^) (Ref. (s^. The molecules crys- 
tallize in a tetragonal structure with lattice parameters 
a = b= 17.319 A and c = 12.388 A. The ground state of 
the molecule has a total electron spin = 10 and, for the 
temperature range of interest in the present work (T < 2 
K), we may describe the electron spin of the cluster by 
means of the effective spin Hamiltonian: 

n = -DSl - BSt + E{Sl - sD ~ C{SX + St) + 

+/iBB.g.S. (1) 

Commonly adopted parameter values are D = 0.548 K, 
B = 1.17 mK and C = 22 /iK as obtained by neu- 
tron scattering dataj^ and for the g tensor the val- 
ues = 1.93 and g± = 1.96 from high-frequency 
EPRi^i^^i^ The uniaxial anisotropy terms —DSl and 
—BSj can be attributed to the single-ion anisotropy of 
the Mn'^+ ions,'^^ which is due to the crystal field effects 
resulting in the Jahn- Teller distortions of the coordina- 
tion octahedra, where the elongation axes are approxi- 
mately parallel to the c-axis of the crystal. Considering 
only the diagonal terms, the energy levels scheme would 
be a series of doublets of degenerate states, | ± m), sep- 
arated by a barrier with a total height DS^ + BS'^ ~ 
66.6 K [Fig. [Hb)]. The transverse anisotropy terms, 
E{Sl - S^) - C{Sl + Si), lift the degeneracy of the 
[ ± m) states and allow quantum tunneling of the giant 



spin through the anisotropy barrier. We call the ma- 
trix element for the tunneling of the giant spin through 
the m-th doublet, and 2 Am the corresponding tunneling 
splitting. The C{S'^ + Si) term arises from the four- 
fold S'4 point symmetry of the molecule, but there is now 
solid experimental evidence^i^^ for the prediction^ that 
a disorder in the acetic acid of crystallization is present 
and gives rise to six different isomers of Mni2 cluster, 
four of which have symmetry lower than tetragonal and 
therefore have nonzero rhombic term E{Sl — Sy). EPR 
experiments give an upper bound E < 14 mKi^ For the 
purpose of NMR experiments, such isomerism may cause 
slight variations in the local hyperfine couplings, caus- 
ing extra broadening in the ^^Mn resonance lines. Very 
recently, a new family of Mni2 clusters has been syn- 
thesized, which does not suffer from the solvent disorder 
mentioned above, and yields indeed more sharply defined 
^^Mn NMR spectra.^" 

When adding spin-phonon interactions, ^^•"'^ the pos- 
sible transitions between the energy levels of ([T|) are 
sketched in Fig. [Tl[b). We distinguish between intrawell 
spin-phonon excitations, where the spin state remains 
inside the same energy potential well, and the interwell 
transitions, which involve spin reversal by quantum tun- 
neling through the barrier, allowed by the terms in ^ 
that do not commute with Sz. Thermally-assisted tun- 
neling involves both these types of transitions. 

The above discussion refers to the majority of the 
molecules in a real sample, but for our experiments the 
crucial feature of Mni2-ac is the presence of fast-relaxing 
molecules (FRMs);^ i.e. clusters characterized by a 
lower anisotropy barrier and a much faster relaxation 
rate, as observed for instance by ac-susceptibility^'* and 
magnetization measurements.^- It has been recognized 
that such FRMs originate from Jahn- Teller isomerism^ 
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i.e. the presence in the molecule of one or two Mn^+ 
sites where the elongated Jahn- Teller axis points in a di- 
rection roughly perpendicular instead of parallel to the 
crystalline c-axis. This results in the reduction of the 
anisotropy barrier to 35 or 15 K in the case of one or two 
flipped Jahn- Teller axes, respectively,— and presumably 
in an increased strength of the non-diagonal terms in the 
spin Hamiltonian as well. Furthermore, the anisotropy 
axis z of the whole molecule no longer coincides with 
the crystallographic c-axis, but deviates e.g. by ~ 10° 
in the molecules with 35 K barrier4^ The Jahn- Teller 
isomerism is very different from the above-mentioned ef- 
fect of disorder in solvent molecules, and produces much 
more important effects for the present study. As will 
be argued below, the presence of the FRMs is essential 
for the interpretation of our results and, to some extent, 
may be regarded as a fortunate feature for this specific 
experiment. 

The sample used in the experiment consisted of about 
60 mg of polycrystalline Mni2-ac, with typical crystallite 
volume ~ 0.1 mm'^. The crystallites were used as-grown 
(i.e., not crushed), mixed with Stycast 1266 epoxy, in- 
serted in a 6 mm capsule and allowed to set for 24 
hours in the room temperature bore of a 9.4 T super- 
conducting magnet. With this procedure, the magnetic 
easy axis of the molecules (which coincides with the long 
axis of the needle-like crystallites) ends up being aligned 
along the field within a few degrees. In addition, we shall 
report NMR spectra taken on a small single crystal (mass 
~ 1 mg). 



B. Low-temperature pulse NMR setup 

Our experimental setup is based on a Leiden Cryogen- 
ics MNK126-400ROF dilution refrigerator, fitted with a 
plastic mixing chamber that allows the sample to be ther- 
malized directly by the ^He flow. A scheme of the low- 
temperature part of the refrigerator is shown in Fig. [21 
together with the NMR circuitry. The mixing chamber 
consists of two concentric tubes, obtained by rolling a 
Kapton foil coated with Stycast 1266 epoxy. The tops 
of each tube are glued into concentric Araldite pots: the 
inner pot receives the downwards flow of condensed '^He 
and, a few millimeters below the inlet, the phase separa- 
tion between the pure '^He phase and the dilute '^He/'^He 
phase takes place. The circulation of '^He is then forced 
downwards along the inner Kapton tube, which has open- 
ings at the bottom side to allow the return of the "^He 
stream through the thin space in between the tubes. 
Both the bottom of the Kapton tail and the outer pot are 
closed by conical Araldite plugs smeared with Apiezon N 
grease. 

A 2-turns copper coil is wound around the capsule con- 
taining the sample, mounted on top of the lower conical 
plug and inserted in the '^He/^He mixture at the bot- 
tom of the mixing chamber tail, which coincides with 
the center of a 9 T superconducting magnet. The coil 



is then connected by a thin brass coaxial cable (length 
K, 0.5 m) to two tunable cylindrical teflon capacitors, 
mounted at the still (see Fig. [2]). At the frequency 
where the cable connecting capacitors and coil is pre- 
cisely one wavelength, the circuit is equivalent to a stan- 
dard lumped LC-resonator. However, since the A-cable is 
a low-conductivity coax for low-T applications, the qual- 
ity factor of the resonator (which includes the cable) is 
drastically reduced. Although this affects the sensitivity 
of the circuit, it also broadens the accessible frequency 
range without the need to retune the capacitors. Cutting 
the cable for one wavelength at ~ 280 MHz, the circuit 
is usable between (at least) 220 and 320 MHz. As for the 
room-temperature NMR electronics, details can be found 
in Ref. Si. 

The temperature inside the mixing chamber is mon- 
itored by two simultaneously calibrated Speer carbon 
thermometers, one in the outer top Araldite pot, and 
the other at the bottom of the Kapton tail, next to the 
sample. At steady state and in the absence of NMR 
pulses, the temperature along the mixing chamber is uni- 
form within < 0.5 mK. The effect of applying high-power 
(~ 100 W) NMR pulses is shown in Fig. [SJa) and (b). 
A sudden increase in the measured temperature is seen 
both at the bottom and the top thermometer, and can 
be attributed to the short electromagnetic pulse. The 
temperature at the lower thermometer, i.e. next to the 
sample and the NMR coil, quickly recovers its unper- 
turbed value, whereas the upper thermometer begins to 
sense the "heat wave" carried by the ^He stream with a 
delay of about 3 minutes. This has the important con- 
sequence that we can use the upper thermometer to dis- 
tinguish the effect of sudden electromagnetic radiation 
bursts from the simple heating of the '^He/^'He mixture, 
as will be shown in ijlVBI below. 

The sample temperature is regulated by applying cur- 
rent to a manganin wire, anti-inductivcly wound around a 
copper joint just above the ^He inlet in the mixing cham- 
ber. In this way we can heat the incoming '^He stream 
and uniformly increase the mixing chamber temperature. 

For the ^He circulation we employ an oil-free pumping 
system, consisting of a 500 ni^/h Roots booster pump, 
backed by two 10 m'^/h dry scroll pumps. The system 
reaches a base temperature of 9 mK, and the practical 
operating temperature while applying r/-pulses is as low 
as 15 - 20 mK. 



C. Measurements and data analysis 

The ^'"'Mn nuclear precession was detected by the spin- 
echo technique. A typical pulse sequence includes a first 
7r/2-pulse with duration 12 — 12 /iS, a waiting interval 
of 45 ytis, and a 24 /is 7r-pulse for refocusing. Given the 
heating effects shown in Fig. [2l a waiting time of 600 s 
between subsequent pulse trains easily allows to keep the 
operating temperature around 15 — 20 mK. Moreover, at 
such low temperature the signal intensity is so high that 
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1 . Coaxial cable connected to 
the NMR spectrometer and 
pulse generator 

2. ^He distillator ("still") with cold plate. 

3. Rotating shafts for the tunable 
capacitors, accessible from the top 
of the refrigerator. 

4. NMR matching capacitor 

5. NMRtuning capacitor. 

6. Upper heat exchanger 

7. 80 mKpot. 

8. ?^-cable, cut to one wavelength 
at 280 MHz. 

9. Lower heat exchanger 

10. Copper joint with heating wire 

11. Araldite mixing chamber 

12. Ur. 

13. Pure ^He phase. 

14. Double-wall Kapton tail. 

15. Forced downwards flow of ^He 
in the dilute phase. 

16. Sample with NMR coil. 

17. Openings in the inner Kapton tube 
to allow the return of the ^He flow. 

18. Lower thermometer. 

19. Conical vacuum plugs. 



FIG. 2: (Color online) Sketch of the low-temperature part of the dilution refrigerator, showing the components of the NMR 
circuitry, the special plastic mixing chamber and the position of the thermometers. Graph panels: temperatures recorded at 
the (a) upper and (b) lower mixing chamber thermometers, having applied a spin-echo NMR pulse sequence at time t = 0. 



we could obtain an excellent signal-to-noise ratio with- 
out need of averaging, so that a typical measurement se- 
quence took less than 12 hours. Above 100 mK it proved 
convenient to take a few averages, but there the heating 
due to the r/-pulses became negligible, and the waiting 
time could be reduced to ~ 100 s. 



The longitudinal spin relaxation (LSR) was stud- 
ied by measuring the recovery of the longitudinal nu- 
clear magnetization after an inversion pulse. We pre- 
ferred this technique to the more widely used saturation 
recovery^iiiSii^ because it avoids the heating effects of 
the saturation pulse train, but we checked at intermedi- 
ate temperatures that the two methods indeed lead to 
the same value of LSR rate. An example of echo sig- 
nals obtained as a function of the waiting time after 
the inversion pulse is shown in Fig. ^a). By integrat- 
ing the echo intensity we obtain the time-dependence of 
the nuclear magnetization, M{t), as shown in Fig. ^h). 
For the ease of comparison between different curves, we 
renormalize the vertical scale such that M{0)/M{oo) — 
-I and M{t > Ti)/M(oo) = 1, even though usually 
|M(0)| < |M(oo)|, as could be deduced from Fig. E^a). 
This is just an artifact that occurs when the NMR line is 
much broader than the spectrum of the inversion pulse, 
and does not mean that the length of the 7r-pulse is in- 
correct. Since the ^^Mn nuclei have spin / = 5/2, we 



fitted the recovery of the nuclear magnetization with 
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where W is the longitudinal spin relaxation rate. Note 
that, in the simple case of a spin 1/2, W is related to 
the relaxation time Ti by 2W = T^^. The above multi- 
exponential expression and its numerical coefficients are 
derived under the assumption that the 1 = 5/2 multiplet 
is split by quadrupolar interactions, and it is possible 
to resolve the central transition within that multiplet. 
While earlier work indicated that all three manganese 
NMR lines are quadrupolar-split,^*^ more recent exper- 
iments on single crystal samples have questioned that 
conclusion, ^^'^^ and thereby the applicability of Eq. ^ 
to the present experiments. Even if other sources of line 
broadening hinder the visibility of the quadrupolar con- 
tribution, the condition for the absence of quadrupolar 
splitting is an exactly cubic environment for the nuclear 
site, which is not satisfied here. For this reason, and for 
the ease of comparison with ou r^'^'^^ and other groups' 
earlier results,— i^*^ we choose to retain Eq. ([2]) for the 
analysis of the inversion recovery data. 

The transverse spin relaxation (TSR) rate was 
obtained by measuring the decay of echo intensity upon 
increasing the waiting time r between the 7r/2- and the 
TT-pulses. The decay of transverse magnetization M±(t) 
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FIG. 3: (Color online) (a) An example of "real time" echo 
signals recorded during an inversion recovery, i.e. measuring 
the echo intensity at increasing delays after an inversion pulse. 
In particular, these are single-shot (no averaging) raw data 
taken at B = and T = 20 mK in the Mn<^' site, (b) The 
(normalized) integral of the echoes (open dots) is fitted to Eq. 
dS]) (solid line) to yield the LSR rate W. 



can be fitted by a single exponential 



M±{2t) 
Mi(0) 



exp 



2t 



(3) 



except at the lowest temperatures (T < 0.2 K), where 
also a gaussian component T^q needs to be included: 



M^{2t) 
Af±(0) 



exp 



2t 



exp 



{2rl 



(4) 



As regards the experiments to determine the nuclear 
spin temperature, the measurements were performed by 
monitoring the echo intensity at regular intervals while 
changing the temperature Tbath of the '^He/'^He bath in 
which the sample is immersed. Recalling that the nuclear 
magnetization is related to the nuclear spin temperature 



Tnuci by the Curie law: 



M(T„uci) = N^lo 



(5) 



and assuming that Tbath = Tnuci at a certain tempera- 
ture To (e.g. 0.8 K), we can define a calibration factor K 
such that Af(To) = if/T'nuci(7o) and use that definition 
to derive the time evolution of the nuclear spin temper- 
ature as Tnuci(0 = K/M{t) while the bath temperature 
is changed. 

Due to the strong magnetic hysteresis of Mni2-ac, it 
is important to specify the magnetization state of the 
sample since, as will be shown below, this parameter can 
influence the observed nuclear spin dynamics. There- 
fore we carried out experiments under both zero-field 
cooled (ZFC) and field-cooled (FC) conditions, which 
correspond to zero and saturated magnetization along 
the easy axis, respectively. Heating the sample up to 
T « 4 K is sufficient to wash out any memory of the 
previous magnetic state. When the sample is already at 
T ^ 1 K, the field-cooling procedure can be replaced by 
the application of a longitudinal field large enough to de- 
stroy the anisotropy barrier, e.g. -B2 = 8 T. importantly, 
the shift of the ^^Mn NMR frequency with external field 
depends on the magnetization state of the sample i^°i^^ 
in a ZFC sample each resonance line splits in two, one 
line moving to + 7n^2 and the other to ujq — 7n-Bz- 
Conversely, in a FC sample only one line is observed, 
shifting to higher or lower frequency depending on the 
direction of Bz relative to the magnetization direction. 
Therefore, by measuring the intensity of the shifted lines 
in a moderate longitudinal field, typically ^ 0.5 T, we 
can check the magnetization of the sample as seen by the 
nuclei that contribute to the NMR signal. 



III. NUCLEAR SPIN DYNAMICS 
A. NMR spectra 

The basic feature of the ^^Mn NMR spectra in Mni2- 
ac is the presence of three well-separated lines, that can 
be ascribed to three crystallographically inequivalent Mn 
sites in the molecule. The Mn^^^ line, centered around 
vi K, 230 MHz, originates from the nuclei that belong 
to the central core of Mn^+ ions, whereas the Mn*^^) and 
Mn^'^) lines, centered at 1^2 ~ 280 and « 365 MHz, re- 
spectively, have been assigned to the nuclei in the outer 
crown of Mn'^^ ions.*^'^" In Fig. [3]we show the Mn^^^ and 
Mn^^) spectra at T = 20 mK, both in the oriented pow- 
der and in the single crystal, in a FC sample. Note that, 
whereas single-crystal spectra of Mni2-ac have been re- 
cently published)^ the present spectra are the only ones 
measured at subkelvin temperatures so far. As argued 
already in Ref. HO, the single-crystal spectra indicate 
that the width of the Mn*^^) line may not originate from 
a small quadrupolar splitting. Instead, at least two in- 
equivalent Mn*+ sites may exist, supporting the growing 
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FIG. 4; (Color online) ^^^Mn NMR spectra of the (a) Mn'^' 
and (b) Mn^^^ lines in Mni2-ac, at T = 20 mK. Open circles: 
oriented powder. Solid squares: single crystal. The Mn sites 
corresponding to each line are shown in the central drawing 
of the molecular structure. All the spectra are measured in a 
field-cooled sample. 



amount of evidence about the lack of symmetry of the 
Mni2-ac compound. 

We also note that the highest peak in the Mn^^' line at 
T = 20 mK is found at a frequency V2 ~ 287 MHz about 8 
MHz higher than most of the previously reported spectra 
at T > 1 K^iiSi^ with the exception of Ref. [H, whereas 

line is consistent with all the 



the position of the Mn^^^ 
previous reports. 



B. Longitudinal spin relaxation in zero field 

The LSR rate as a function of temperature for the 
Mn^^) line, in zero field and zero-field cooled (ZFC) sam- 
ple, is shown in Fig. [S] The most prominent feature in 
these data is a sharp crossover at T ~ 0.8 K between 
a roughly exponential T-dependence and an almost T- 
independent plateau. We have previously attributed the 
T-independent nuclear relaxation to the effect of tunnel- 
ing fluctuations within the ground doublet of the cluster 
spinS)^ and we shall dedicate most of the present pa- 
per to discuss our further results supporting this state- 
ment. Here we shall also argue that, even in the high- 
temperature regime, thermally assisted quantum tunnel- 
ing plays an essential role, and the experimental results 
cannot be understood simply in terms of LSR driven by 
intrawell electronic transitions.'*^ It should be noted that 
the crossover from thermally activated to ground-state 
tunneling has also been observed by analyzing the T- 
dependence of the steps in the magnetization hysteresis 
loops i^i^ The important advantage of our NMR mea- 
surements is that the nuclear dynamics is sensitive to 
fluctuations of the cluster electron spins without even re- 
quiring a change in the macroscopic magnetization of the 
sample. Clearly, no macroscopic probe (except perhaps 
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FIG. 5: (Color online) Temperature-dependence of the nu- 
clear spin-lattice relaxation rate W of the Mn'^' line, in zero 
external field and ZFC sample. The inset shows some exam- 
ples of recovery of the nuclear magnetization after a time t 
from an inversion pulse, at the indicated temperatures. These 
curves have been fitted to Eq. ([2} to extract W . 



an extremely sensitive magnetic noise detector) would 
be able to detect the presence of tunneling fluctuations 
in a zero- field cooled sample in zero external field, since 
the total magnetization is zero and remains so. Below 
T ~ 1.5 K the steps in the hysteresis loops of Mni2-ac 
can be observed only at relatively high values of external 
fieldf^i^ which means that the spin Hamiltonian under 
those conditions is radically different from the zero-field 
case. Therefore, that both our data and the previous 
magnetization measurements show a crossover around 
T ~ 0.8 K should be considered as a coincidence. 

The roughly T-independent plateau in the LSR rate 
below T ~ 0.8 K is characterized by a value of ~ 0.03 
s~* which is surprisingly high, which at first sight may 
appear like an argument against the interpretation in 
terms of tunneling fluctuations of the electron spin. Ex- 
perimentally it is indeed well knownS that the relaxation 
of the magnetization in Mni2-ac in zero field may take 
years at low T, which means that the tunneling events 
are in fact extremely rare. Based on this, we are forced to 
assume that tunneling takes place only in a small minor- 
ity of the clusters, and that some additional mechanism 
takes care of the relaxation of the nuclei in molecules that 
do not tunnel. This is a very realistic assumption, since 
all samples of Mni2-ac are reported to contain a fraction 
of FRMsf^i^ as mentioned in Sect. Ill Al Moreover, since 
we are also able to monitor the sample magnetization, we 
verified that e.g. a FC sample maintains indeed its satu- 
ration magnetization for several weeks while nuclear re- 
laxation experiments are being performed (at zero field). 



This confirms that any relevant tunnehng dynamics must 
originate from a small minority of molecules. On the 
other hand, it also means that the observed NMR signal 
comes mainly from nuclei belonging to frozen molecules, 
thus there must be some way for the fluctuations in FRMs 
to influence the nuclear dynamics in the majority of slow 
molecules as well. One possibility is to ascribe it to the 
fluctuating dipolar field produced by a tunneling FRM at 
the nuclear sites of neighboring frozen molecules. In that 
case we may give an estimate of W using an expression 
of the form: 



4 dip-j^ 



Tt 



-'N'T 



(6) 



Trj7^ 3> 10^ s^^. Such a high rate is of course 
completely unrealistic. We must therefore consider the 
effect of a tunneling molecule on the nuclei that belong to 
the molecule itself, and look for some additional mecha- 
nism that links nuclei in FRMs with equivalent nuclei in 
frozen clusters. It is natural to seek the origin of such a 
mechanism in the intercluster nuclear spin diffusion, and 
in the next section we shall provide strong experimental 
evidences to support this interpretation. 
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where b^ip is the perpendicular component of the fluctu- 
ating dipolar field produced by a tunneling molecule on 
its neighbors and t^^ is the tunneling rate. The highest 
value that &dip may take is 3 mT in the case of near- FIG. 6: 
est neighbors, which leads to the condition W ~ 0.03 rate T,' 



0.01 



0.1 



T (K) 



(Color online) Temperature-dependence of the TSR 
^ (squares) rates for a ZFC sample in zero field and 



u — 231 MHz. The solid line in the T-independent regime 
is a guide for the eye. Inset: normalized decay of transverse 
nuclear magnetization, M{t)/M{0), for ZFC (full squares) 
and FC (open squares) sample, at T = 20 mK. The solid lines 
are fits to Eq. Q, yielding the ratio r2c^(FC)/r2c^(ZFC) = 
1.35 ~ \/2. The sketches in the inset represent pictorially the 
fact that intercluster spin diffusion is possible in a FC sample 
since all the nuclei have the same Larmor frequency, contrary 
to the case of a ZFC sample. 



C. Transverse spin relaxation 

The T-dependence of the TSR rate T^^{T) is shown 
in Fig. [SI One may observe that below 0.8 K the TSR, 
just like the LSR, saturates to a nearly T-independent 
plateau. In particular, r2"^(T < 0.8 K) w 100 s^^, 
which is a factor ~ 3000 larger than the low-T limit of the 
LSR rate W. The values plotted in Fig. [Slare all obtained 
by fitting the decay of the transverse magnetization with 
Eq. ([3]), i.e. with a single exponential. While this is very 
accurate at high T, we found that for T < 0.2 K a better 
fit is obtained by including a Gaussian component, as 
in Eq. (|4]). In any case, the single-exponential fit does 
capture the relevant value for at all temperatures. 

A point of great interest is the measurement of the 
TSR at T = 20 mK in a FC and a ZFC sample, as shown 
in the inset of Fig. [6l The decay of the transverse magne- 
tization is best fitted by Eq. ([U, whereby the Gaussian 
component, , is separated from the Lorentzian one, 
T^i^ . From the Gaussian component of the decay we can 
extract directly the effect of the nuclear dipole-dipole in- 
teraction, whereas the other mechanisms of dephasing 
(e.g. random changes in the local field due to tunneling 
molecules) contribute mainly to the Lorentzian part. The 
fit yields T-(}{FC) = 104±3 s'^ and T-(}{ZFC) = 77 ±3 
s^^. These results can be understood by assuming that, 
at very low T, the main source of TSR is the dipole- 
dipole coupling of like nuclei in neighboring molecules. 



Then we can estimate from the Van Vleck formula 
for the second moment M2 = (Aw^) of the absorption 
line in dipolarly-coupled spinsi^ 

i>j V 

yielding T^^ — 131 s~^ if we take for Vij the distance be- 
tween centers of neighboring molecules. The estimated 
Tj"^ would obviously be much larger if one would consider 
the coupling between nuclei within the same cluster. As 
we argued when discussing the ^^Mn spectra, it is possi- 
ble that the cluster symmetry is low enough to prevent 
intracluster nuclear spin fiip-flops. This may explain why 
Eq. ([7]) yields the right order of magnitude when only 
coupling between nuclei in neighboring molecules is con- 
sidered. An alternative argument is that, given the small 
number (4 at best) of like ^^Mn spins within one clus- 
ter, the dipolar coupling between them does not yield a 
genuine decay of the transverse magnetization for the en- 
tire sample. The macroscopic T2 decay measured in the 
experiment reflects therefore the slower, but global, in- 
tercluster spin diffusion rate. A similar observation was 
recently made also in a different molecular compound, 
AI50C120H180 (Ref. M). 

We also note that, in the case of a ZFC sample, the 
sum in Eq. ([7]) should be restricted to only half of the 
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neighboring molecules, since on average half of the spins 
have resonance frequency +wn smd the other half —uj-^^, 
and no flip-flops can occur between nuclei experiencing 
opposite hyperfine fields. This is equivalent to diluting 
the sample by a factor 2, which reduces the expected 
in ZFC sample by a factor \/2. Indeed, we find in 
the experiment T2Q^(FC)/r2G^(ZFC) = 1.35 ~ ^/2 which, 
together with the good quantitative agreement with the 
prediction of Eq. ([7]), constitutes solid evidence for the 
presence of intercluster nuclear spin diffusion. This is 
precisely the mechanism required to explain why the tun- 
neling in a minority of FRMs can relax the whole nuclear 
spin system. The need for intercluster nuclear spin dif- 
fusion could already have been postulated by analyzing 
the LSR rate, and the magnetization dependence of the 
TSR rate gives an independent confirmation. 

For comparison, in a recent study of the ^^Fe NMR 
in Fcs, Back et al.^^ attributed the observed TSR rate 
to the dipolar interaction between ^^Fe and nuclei. 
They analyzed their data with the expression Tj^^ ~ 

(M2^"Vl2rc)^/^, where is the proton TSR time due 

(a) 

to their mutual dipolar coupling and M2 is the sec- 
ond moment of the '^^Fe - coupling. However, the 
same model^ predicts the echo intensity to decay as 
M^(i)/M^(0) ~ exp(-2M2^"^t3/3). This function fails 
completely in fitting our echo decays, therefore we do not 
consider the ^^Mn - dipolar coupling as an alternative 
explanation for the TSR we observe. 

Finally we stress that, in our view, the fact that the 
LSR and the TSR are both roughly T-independent be- 
low 0.8 K does not find its origin in the same mechanism. 
Rather, we attribute them to two different mechanisms, 
both T-independcnt: the quantum tunneling of the elec- 
tron spin (for the LSR) and the nuclear spin diffusion 
(for the TSR). 

Having argued that the LSR in Mni2-ac in driven by 
tunneling fluctuations of the FRMs, which are peculiar 
of the acetate compound, it's interesting to note that 
other varieties of Mni2 molecules have meanwhile be- 
come available. In particular the Mni2-ffiuAc^^'§i is a 
truly axially symmetric variety that does not contain any 
FRMs, and could provide an interesting counterexample 
for our results if studied by low-T NMR. The Mni2BrAc 
molecule is also thought to be free of FRMs,— and some 
low-T NMR experiments have been performed on it^I 
that show indeed very different results from what we re- 
port here. However, as we shall argue in mVl a definite 
conclusion on the meaning of NMR experiments at very 
low T should only be drawn when the analysis of the 
nuclear spin thermalization is included. 



D. Field dependence of the longitudinal spin 
relaxation rate 

Further insight in the interplay between the quantum 
tunneling fluctuations and the nuclear spin dynamics is 
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FIG. 7: (Color online) Longitudinal field dependence of the 
LSR rate W in the ZFC (solid dots) and FC (open dots) 
sample at T = 20 mK. The measuring frequencies are i^{Bz) = 
230 + jnBz MHz. The solid line is a Lorentzian fit with 
HWHM AB^ ~ 60 mT. The dotted line through the FC data 
is a guide for the eye. 



provided by the study of the dependence of the LSR on 
a magnetic field Bz applied along the anisotropy axis. It 
is clear from the Hamiltonian ([T]) that, in the absence 
of other perturbations, such a longitudinal field destroys 
the resonance condition for electron spin states on oppo- 
site sides of the barrier and therefore inhibits the quan- 
tum tunneling. In the presence of static dipolar fields, 
-Bdip, by studying the tunneling rate as a function of 
Bz one may in principle obtain information about the 
distribution of longitudinal i?dip, since at a given value 
of Bz there will be a fraction of molecules for which 
_B(jip — —Bz and will therefore be allowed to tunnel just 
by the application of the external bias. 

We show in Fig. [7]the LSR rate W{Bz) at T = 20 mK 
in the ZFC sample, obtained while shifting the measure- 
ment frequency as v{Bz) = v{{)) + ^nBz with v{{)) = 230 
MHz, in order to stay on the center of the NMR line 
that corresponds to the molecules that are aligned ex- 
actly parallel with the applied field. Since for a ZFC 
sample the magnetization is zero, the field dependence 
should be the same when Bz is applied in opposite di- 
rections, as is observed. The data can be fitted by a 
Lorentzian with a half width at half maximum (HWHM) 
/S.Bz — 60 mT: this differs both in shape (Gaussian) 
and in width (Ai?^ ~ 21 mT) from the calculated dipo- 
lar bias distribution in a ZFC sample. An alternative 
experimental estimate, Ai?^ ~ 25 mT, can be found in 
magnetization relaxation experiments,'*^ but only around 
the first level crossing for FRMs (~ 0.39 T) in the FC 
sample. For comparison. Fig. [7] also shows W{Bz) in 
the FC sample: the shape is now distinctly asymmetric, 
with faster relaxation when the external field is opposed 
to the sample magnetization. Interestingly, W{Bz) in 
the FC sample falls off much more slowly on the tails 
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for both positive and negative fields, while the value at 
zero field is less than half that for the ZFC sample. We 
therefore observe that in zero field the recovery of lon- 
gitudinal magnetization in the FC sample is faster than 
in the ZFC, whereas the opposite is true for the decay of 
transverse magnetization (inset Fig. [6|) . 

If the LSR rate W{Bz) is to be interpreted as a signa- 
ture of quantum tunneling, its HWHM is clearly larger 
than expected. Part of the reason may be the fact that 
the width of the Mn^^^ line is already intrinsically larger 
than both ABz and the distribution of dipolar fields cre- 
ated by the molecules. Indeed, the width of the Mn*^^^ 
line, a„ ~ 1.2 MHz, translates into a local field distri- 
bution of width (Tb — 115 mT for ^'"^Mn. The observed 
HWHM does depend, for instance, on the choice of i^{0). 
As soon as Bz 7^ the presence of slightly misaligned 
crystallites in our sample may also contribute to the 
width of the resonance. In any case, all of the mecha- 
nisms mentioned above (distribution of internal dipolar 
fields, width of the NMR line, distribution of crystallite 
orientations in the sample) would yield a T-independent 
linewidth for W{B,). Fig. [8] shows W{B,) in ZFC sam- 
ple at three different temperatures, T — 0.02,0.72,1.13 
K, covering the pure quantum regime, the thermally- 
activated regime, and the crossover temperature. The 
NMR frequency in these datasets is v{Bz) = 22,l + "fNBz- 
The data have been fitted by Lorentzian lines yielding a 
HWHM AB^ = 16,85, 118 increasing with temperature. 
We note immediately that the HWHM at T = 20 mK 
is much smaller than the one obtained from the data in 
Fig. [71 the only difference between the two sets being v{{)) 
and, subsequently, all other measurement frequencies at 
Bz ^ 0. Indeed, we found that already in zero field the 
LSR rate does depend on v, reaching the highest values 
at the center of the line and falling off (up to a factor 5) 
on the sides. This dependence, however, becomes much 
weaker at high temperatures. It is therefore rather diffi- 
cult to make strong statements about the meaning of the 
observed increase in Ai?z with temperature. At any rate, 
however, the field dependencies observed here at low-T 
are much stronger than those previously reported in the 
high-T regime<^i^ Goto et al. also reported W{Bz) for 
the "lower branch" of the Mn^^^ line, viz. for the nuclei 
whose local hyperfine field is opposite to the external field 
(Ref. 24, Fig. 6, closed squares). That situation is equiv- 
alent to our FC data (Fig. [71 open dots) for Bz < 0. At 
large fields an overall increase of W with Bz is observed 
in Ref. [l^, but for < 1 T the LSR rate does decrease, 
in agreement with our results. 

We also noted, both in Fig. [51 and in the FC data 
in Fig. [71 that a small increase in W{Bz) occurs at 
\Bz\ — 0.5 T, which is approximately the field value at 
which the | -I- 9) and | — 10) electron spin states come 
into resonance. This feature is barely observable, but 
nevertheless well reproducible. As a counterexample, in 
another dataset (not shown) we investigated W{Bz) more 
carefully in the FC sample at T = 20 mK for positive val- 
ues of Bz, and found no increase around Bz — 0.5 T, as 
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FIG. 8: (Color online) Longitudinal field dependence of the 
LSR rate in ZFC sample at T = 20 mK (down triangles), 
T = 720 mK (diamonds) and T = 1.13 K (up triangles). The 
measuring frequency in these datasets is = 231 -I- ^nBz 
MHz. The lines are Lorentzian fits yielding HWHM AB^ = 
16, 85, 118 mT, respectively. 

one would expect since the fully populated state, | — 10), 
is pushed far from all other energy levels. A similarly 
small peak in W{Bz) at the first levels crossing has been 
recently observed in Fes as welli^ 



E. Deuterated sample 

The role of the fluctuating hyperfine bias on the in- 
coherent tunneling dynamics of SMMs, predicted by 
Prokof'ev and Stamp, ^ has been clearly demonstrated 
by measuring the quantum relaxation of the magnetiza- 
tion in Feg crystals in which the hyperfine couplings had 
been artificially modified by substituting '^^Fe by ^^Fe or 
^H by ^H (Ref. [isl ). For instance, the time necessary 
to relax 1% of the saturation magnetization below 0.2 
K was found to increase from 800 s to 4000 s by sub- 
stituting protons by deuterium, whereas it decreased to 
300 s in the ^^Fe enriched sample. More recently, Evan- 
gelisti et al)^ showed that the ^^Fe isotopic enrichment 
of Fcg causes the magnetic specific heat to approach its 
equilibrium value within accessible timescales (~ 100 s). 

Since in Mni2-ac the only possible isotope substitu- 
tion is ^H ^^H, we performed a short set of measure- 
ments on a deuterated sample. The sample consists of 
much smaller crystallites than the "natural" ones used 
in all other experiments reported here. Although a field- 
alignment was attempted following the same procedure 
as described in i jllCl the orientation of the deuterated 
sample turned out to have remained almost completely 
random, probably due to the too small shape anisotropy 
of the crystallites. We therefore report only experiments 
in zero external field, where the orientation is in principle 
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FIG. 9; (Color online) Comparison between (a) the nuclear inversion recoveries and (b) the decays of transverse magnetization 
in the "natural" Mni2-ac (circles) and in the deuterated sample (squares), at T = 20 mK in zero field and ZFC sample, for the 
Mn'^' site. The solid lines in (a) are fits to Eq. (gj. 



irrelevant. 

The results are shown in Fig. [HI the ^^Mn LSR rate 
at T = 20 mK in zero field and ZFC sample is indeed re- 
duced to Wdeut — 0.0035 s~^, i.e. 6.5 times lower than in 
the "natural" sample. This factor is the same as the re- 
duction of the electron spin relaxation rate seen in deuter- 
ated Fes (Ref. M), and it coincides with the ratio of the 
gyromagnetic ratios of and ^H. This finding unequiv- 
ocally proves that the proton spins are very effective in 
provoking the tunneling events via the Prokof'ev-Stamp 
mechanism, and confirms that the LSR rate of the ^^Mn 
nuclei is a direct probe of the electron spin tunneling rate. 

As regards the TSR, the result is quite intriguing: slow 
but rather ample oscillations are superimposed to the de- 
cay of transverse magnetization, and the overall decay 
rate appears slower than in the natural sample. This be- 
havior is reminiscent of the change in TSR rate upon ap- 
plication of a small longitudinal magnetic field in the nat- 
ural sample. The latter has a rather complicated physical 
origin and is still under investigation. 

F. Comparison with a Mn"'"'" site 

Some rather interesting results emerge from the anal- 
ysis of extra measurements performed on the NMR line 
of the Mn^^^ site, i.e. a Mn'^+ ion. Fig. [TUl shows a com- 
parison between the recovery of the longitudinal magne- 
tization and the decay of the transverse magnetization in 
Mn(i) and Mn^^) sites, at T = 20 mK in the FC sample 
and zero external field, at a frequency j/*^^) = 283.7 MHz. 
The TSR is very similar in both sites, although a closer 
inspection evidences that the Gaussian nature of the de- 
cay is less pronounced in the Mn*^^) sites, which leads to 
T-(^ = 83 s-i instead of the T'^} = 104 s'^ found in 
Mn*^^-'. More importantly, the LSR is three times slower 
in the Mn^^^ site, as seen in Fig. [TUT b). This is oppo- 
site to the high-T regime, where the Mn'^+ sites were 
foundSli^ to have much faster relaxation. Furthermore, 



the field dependence of the LSR rate appears sharper in 
the Mn^^^ site, as shown in Fig. [TT] The asymmetry 
in W{Bz) for a FC sample is still present, but less evi- 
dent than in the Mn^^^ site due to the more pronounced 
decrease of W already for small applied fields. 

The similarity between the TSR rates in the Mn^^^ and 
the Mn*^^) sites is indeed expected if T2 is determined 
by intercluster nuclear spin diffusion. Conversely, the 
difference in LSR is more difficult to understand if one 
assumes that the process that induces longitudinal spin 
relaxation is the tunneling of the molecular spin. How- 
ever, one clear difference between Mn^^^ and Mn^^-* is the 
width of the NMR line, much larger in Since the 

integrated intensity of both lines is identical, the Mn'^' 
has an accordingly lower maximum intensity. We have 
verified for both sites that the LSR rate is the fastest 
when measuring at the highest intensity along each line. 
Thus, the factor 3 slower LSR in Mn*^^^ could simply be 
another manifestation of the apparent dependence of the 
measured W on the NMR intensity along each line. We 
point out, however, that the measured LSR rate is inde- 
pendent of the 7r/2 pulse length, which determines the 
spectral width of the pulse and thereby the fraction of 
spins being manipulated and observed. This means that 
the difference in W for the two sites cannot be simply 
attributed to a difference in the number of spins excited 
during a pulse of given length but that other (more com- 
plex) factors must play a role. 



IV. THERMALIZATION OF THE NUCLEAR 
SPINS 

Having demonstrated that the ^^Mn longitudinal spin 
relaxation below 0.8 K is driven by T-independent quan- 
tum tunneling fluctuations, a natural question to ask is 
whether or not the nuclear spins are in thermal contact 
with the lattice at these low temperatures. Let us re- 
call that any direct coupling between phonons and nu- 
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FIG. 10; (Color online) Comparison between (a) the recovery of longitudinal magnetization and (b) the decay of transverse 
magnetization in Mn*-^^ (circles) and Mn^^' (diamonds) sites, at T = 20 mK in FC sample and zero external field. The solid 
(Mn'^') and dashed (Mn'^') lines are fits to Eq.© in panel (a) and Eq.Q in panel (b). 



clear spins is expected to be exceedingly weak, due to 
the very small density of phonons at the nuclear Larmor 
frequencyi^ Relaxation through electric quadrupole ef- 
fects, if present, would show a temperature dependence 
cx (T/9d) for direct process or cx (T/9d)^ for Raman 
process (0d is the Debye temperature), which is not con- 
sistent with our observations. Therefore the thermaliza- 
tion of the nuclei will have to take place via the electron 
spin - lattice channel. Since in the quantum regime the 
only electron spin fluctuations are due to tunneling, the 
question whether the nuclear spins will still be in equi- 
librium with the lattice temperature is of the utmost im- 
portance. 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

2-160 

E 140 

"^120 
<u 

i: 100 




(0 80 




3 6 9 12 15 
Time (h) 



2 3 4 

Time (h) 




-0.2 0.0 



FIG. 11: (Color online) Longitudinal field dependencies of 
the LSR rates in Mn'^' (circles) and Mn'^' (diamonds) sites, 
normalized at the zero-field value. The data are taken at 
T = 20 mK in FC sample with central measuring frequencies 
v'-^\0) = 230 MHz and i/'^'(0) = 283.7 MHz. 



FIG. 12: (Color online) Comparison between bath temper- 
ature Tbath (solid lines) and nuclear spin temperature Tnuci 
(circles), while cooling down the system (main panel) and 
while applying step-like heat loads (inset). The waiting time 
between NMR pulses was 60 s in the main panel and 180 s in 
the inset. Both datasets are at zero field in ZFC sample. 



A. Time evolution of the nuclear spin temperature 

We have addressed this problem by cooling down the 
refrigerator from 800 to 20 mK while monitoring simul- 
taneously the temperature Tbath of the '^He/^He bath in 
the mixing chamber (just next to the sample) and the 
NMR signal intensity of the Mn^^^ line, in zero external 
field and on a ZFC sample. The signal intensity was mea- 
sured by spin echo with repetition time irop = 60 s. The 
nuclear spin temperature^ Tnuci is obtained as described 
in OTCl and plotted in Fig. [HI together with Tbath- We 
find that the nuclear spin temperature strictly follows 
the bath temperature, with small deviations starting only 
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TABLE I: Experimental conditions and relaxation rates for 
the nuclear spin temperature experiments in Fig. [TJ] 



phenomenological function: 



Panel 


Mn 




h 


Q 


^rcp 


Ti 


Tth 




site 


(T) 


(/imol/s) 


(mW) 


(s) 


(s) 


(min) 


a 


1 





330 


0.63 


120 


41.3 


58 ±5 


b 


2 





330 


0.63 


120 


122 


83± 13 


c 


1 





430 


0.78 


120 


41.3 


37 ±3 


d 


2 


0.2 


330 


0.63 


300 


355 


92 ±33 



below ~ 200 mK. This result is crucial but rather para- 
doxical, and we shall discuss its implications in detail in 
PVDI Experimentally, however, it certifies the effective- 
ness of our cryogenic design in achieving the best possible 
thermalization of the sample, since the nuclear spins are 
the last link in the chain going from the ^He/''He bath 
via the phonons in the sample to the electron spins and 
finally to the nuclei. 

The lowest spin temperature that can be measured ap- 



pears to depend on the pulse repetition time t 



rep ■ 



To 

measure Tnuci with the pulse NMR method we need a 
tt/2 pulse to create a transverse nuclear magnetization, 
and after a time T2 the spins are effectively at infinite 
T so enough time must elapse before taking the next 
Tnuci measurement. For the data in Fig. [T^l irep — 60 
s was barely longer than the observed time for inversion 
recovery [see Fig. EJb)], and the lowest observed spin 
temperature is T^JJ ~ 80 mK. This improved when us- 
ing longer waiting times between pulses, e.g. T^JI ~ 35 
mK with ircp = 180 s, as shown in the inset of Fig. [T^ 
However, no matter how long the waiting time, we never 
observed a Tnuci lower than ~ 30 mK. 

Next we study the time constant Tth for the thermal- 
ization of the nuclear spin system with the helium bath, 
by applying step-like heat loads and following the time 
evolution of Tnuci- In particular, we are interested in the 
relationship between Tth, the LSR time Ti = 1/2W as ob- 
tained from the inversion recovery technique, and the '^He 
circulation rate n, which is proportional to the refriger- 
ator's cooling power, Q. Ti is easily tuned by measuring 
at different longitudinal fields and Mn sites, while fi is 
changed by applying extra heat to the refrigerator still. 
Since also the NMR signal intensity changes under differ- 
ent fields and Mn sites, we must redefine every time the 
conversion factor K between signal intensity and Tnuci- 
In the following we choose K such that the asymptotic 
value of Tnuci for t — > 00 matches the measured Tbath 
at the end of the heat step. This implies the assumption 
that the measuring pulses do not saturate, i.e. "heat up" , 
the nuclear spins, and requires Tiop > Ti. Fig. [T^ shows 
four examples of the time evolution of Tnuci under the 
application of a heat load for ~ 2 hours, in Mn*^^^ and 
Mn^^) sites, with or without an applied field, and with 
an increased ^He flow rate. We fitted the data to the 



[Tnuci (00) - Tn, 



Tnucl(^) — Tnucl(O) -f 
Tth 



(8) 



1 — exp — 



where Tnuci(oo) is set by definition equal to Tbath at the 
end of the step, Tnuci (0) follows automatically from the 
above constraint, and to is the time at which the heat 
pulse is started. We find that Tth is always much longer 
than the nuclear LSR time Ti , and that larger Ti corre- 
sponds to larger Tth- However, the dependence of Tth on 
Mn site and applied field is not as strong as for Ti , i.e. Tth 
and Ti are not strictly proportional to each other. Con- 
versely, by changing the ''He flow rate we observe that, 
within the errors, the ratio of heat transfer from the '^He 
stream to the nuclear spins is proportional to n, given 
the same conditions of nuclear site and external field. 

We should stress that, when measuring Ti by inversion 
recovery, we effectively "heat up" only a small fraction of 
the nuclear spins, namely those whose resonance frequen- 
cies are within a range, Sv, proportional to the inverse 
of the duration, t^, of the 7r-pulse. With ~ 20 /is we 
get Siy = l/27rt7r — 8 KHz, which is less than 0.2% of 
the width of the Mn^^^ line. Conversely, by increasing 
the bath temperature we heat up the entire spin sys- 
tem, thereby requiring a much larger heat flow to occur 
between the ''He stream and the nuclear spins. There- 
fore, these results show that the thermal equilibrium be- 
tween nuclear spins and lattice phonons does occur on a 
timescale of the order of Ti as obtained from inversion 
recovery, since the main bottleneck appears to be be- 
tween lattice phonons and ''He stream, as demonstrated 
by the dependence of Tth on h. In a later set of ex- 
periments (not shown here) using a small single crystal 
instead of a large amount of oriented powder, we have in- 
deed observed an even shorter Tth, which indicates that 
Tth should ultimately tend to Ti for small sample size 
and strong thermal contact between lattice phonons and 
helium bath. 



Longitudinal field sweeps and magnetic 
avalanches 



To conclude our study on the nuclear spin thermal- 
ization, we attempted to measure Tnuci in the presence 
of large longitudinal magnetic field sweeps, motivated by 
the fact that much of the experiments on spin tunneling 
in SMMs are based on the measurement of magnetic hys- 
teresis loops. Under those conditions, the electron spins 
are flipped at abnormally large rates, and one may ask 
whether or not the nuclear spins are still able to remain 
in thermal equilibrium. Unfortunately, monitoring Tnuci 
while Bz is being swept means that one should continu- 
ously change the NMR probe frequency, and synchronize 
that change with the field sweep. This being technically 
cumbersome, we could only mecLSure Ximci zero field. 
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FIG. 13: (Color online) Time evolution of the nuclear spin temperature (open symbols) and the bath temperature (dotted 
lines) upon application of a step-like heat load. All data are for a FC sample. The solid lines are fits to Eq. (8), yielding the 
thermal time constants Tth reported in table [H along with the Mn site, external magnetic field B^, LSR rate W, NMR pulse 
repetition time Uep, ^He fiow rate h, and applied heat load Q. Notice in particular the effect of a change in ^He circulation 
rate, panel (c) vs. panel (a). 



at the beginning and at the end of a Bz sweep. The 
results are somewhat inconclusive and shall not be dis- 
cussed here, but more details can be found in section 
4.4.2 of Ref. El. 

We do mention, however, that during the i3z-sweep ex- 
periments we always encountered magnetic avalanches, 
i.e., abrupt reversal of the electronic magnetization of 
the whole sample. This phenomenon has been first re- 
ported already some time agc^^ but is only recently being 
studied in more detail. Importantly, the magnetization 
reversal is expected to be accompanied by the emission 
of electromagnetic radiation;^ which is in fact what we 
observed in our experiments, since we were not equipped 
to measure the electronic magnetization directly on short 
time scales. Fig. [T4l shows the temperature recorded by 
the upper thermometer in the mixing chamber (see Fig. 
(2), while the longitudinal field is being swept at a rate 
dBz/dt = 0.5 T/min. The sweeping field gives a heat 
load that raises the observed temperature to ~ 30 mK, 
but the most striking feature of the data is the sudden 
jump of Tuppor to above 100 mK, whenever the applied 
field reaches \Bz\ ~ 1.9 T and its direction is oppo- 
site to the instantaneous magnetization. We note that 
the timescale for the apparent temperature jump is es- 
sentially identical to what we observe immediately after 
the application of a rf-pulse for NMR measurements, as 



shown in Fig. ^h). In the same figure it is seen that 
a heat pulse applied at the sample location shows its ef- 
fect at the upper thermometer with a delay of about 3 
minutes (due to the "^He drift velocity) in the form of 
a broad temperature "bump" . We therefore conclude 
that the sudden jumps in Tuppcr shown in Fig. IMf b) 
must be of electromagnetic rather than thermal origin, 
and may be attributed to the radiation produced by the 
sudden reversal of the entire electronic magnetization of 
the sample by the magnetic avalanche^SS, The radiation 
bursts reported in Ref. ^68 at a temperature T — 1.8 
K occurred at ~ 1.4 T, which corresponds to 

the third level crossing field for spin tunneling, i.e., the 
value of field at which the resonance between m = ±10 
and TO — =f7 states is obtained. We found instead the 
avalanches at Isi^'^^l ~ 1.9 T, i.e., the fourth level cross- 
ing, TO = ±10 <-> =f6, but our measurements are done at 
T ~ 30 mK. Goto et al^ also reported the observation of 
magnetic avalanches in Mni2-ac, and studied the temper- 
ature dependence of the avalanche field Bi^'^^ . Their find- 
ing that si''^^ increases with temperature was interpreted 
as a sign that the avalanches occur more easily when the 
thermal contact to the bath is weaker. Indeed, whereas 
they would observe avalanches even at fields as low as 
^(av) 2± 0.5 T (the first level crossing) with the sample 
loosely anchored to the mixing chamber of a dilution re- 
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FIG. 14: (Color online) (a) Longitudinal magnetic field and 
(b) temperature of the upper thermometer (see Fig. [2J during 
a field sweep at dB^/dt — 0.5 T/min. The sample was initially 
field-cooled with Bz > 0. The sharp jumps in Tupper occur 
when \Bz\ — 1.9 T, i.e. at the fourth level crossing field, 
and are attributed to the radiation produced by a magnetic 
avalanche. 



frigerator atT — 0.15 K, they never saw avalanches when 
the same sample was placed directly in a liquid helium 
bath at T = 1.4 K. In this sense, our observation of a 
high Bi''"^ ~ 1.9 T confirms once more that our strategy 
for the sample thermalization is very effective. Suzuki et 
al^ found even higher values of B^^^ at subkelvin tem- 
peratures when measuring the local magnetization of a 
small Mni2-ac crystal immersed in liquid '^He. However, 
their observations differ markedly from ours in that they 
found avalanches occurring in a wide range of (not neces- 
sarily resonant) fields, whereas we saw avalanches always 
and only at the fourth level-crossing field. 



V. ANALYSIS OF THE NUCLEAR SPIN 
DYNAMICS AND THEORETICAL 
IMPLICATIONS 

In this section we attempt a quantitative analysis of 
our experimental results, particularly the observed values 
of LSR rate. To this end, we shall apply the Prokof'ev- 
Stamp (PS) theory of the spin bath, which describes the 
dynamics of a "central spin" S (here the giant electronic 
spin of a Mni2-ac cluster) coupled to a bath of envi- 
ronmental (in this case, nuclear) spins. In view of the 
complexity of the model we provide here an introduc- 
tory overview of some essential elements of the PS the- 
ory needed for our analysis, referring the reader to the 
original papers^iii^S for more details. For comparison. 



we also calculate the LSR rate assuming that the electron 
spin tunneling is driven by spin-phonon coupling JSili- We 
anticipate that the result of this effort will be that the ex- 
isting theory is not sufficient to properly describe these 
and other related experimentsJ^i^^ We shall carry out 
the analysis in detail in order to emphasize at every step 
what assumptions are being made, what is their actual 
validity, and why the known theories cannot explaining 
the data. 



The goal of our analysis is to link the electron spin 
tunnehng rate, F = t^^ , to the observed LSR rate, W, 
based on the following assumptions, justified by the ex- 
periments presented in the previous sections: (i) The nu- 
clear relaxation is driven by tunneling fiuctuations in a 
minority of fast-relaxing molecules. We shall assume the 
fraction of FRMs to be 5% of the total4^ The neighboring 
slow molecules can be safely considered as frozen during 
the timescale of interest and serve simply as a "reservoir 
of nuclear polarization" . (ii) The dipole-dipole coupling 
between ^^Mn nuclei in equivalent sites of neighboring 
molecules allows intercluster nuclear spin diffusion, at a 
rate T^^ much faster than the LSR rate, (iii) The nuclear 
spin system is in thermal equilibrium with the phonon 
bath. 



Before we start, it is of interest to point out some 
rather striking peculiarities of the problem at hand. 
First and most importantly, one cannot use any result 
from perturbation theory here, because the nuclear Zee- 
man splittings arise uniquely from hyperfine fields, which 
themselves jump between two different directions each 
time the electron spin of a molecule tunnels, so there is no 
static part of the nuclear Hamiltonian. Perhaps the only 
situation that resembles this is the nuclear quadrupolar 
relaxation in systems with molecular rotations. ^"^ Con- 
versely, in the overwhelming majority of NMR experi- 
ments one has a static external field (produced by an 
actual magnet) and some local fluctuating fields arising 
from the magnetic environment of the nuclei, which can 
be treated as small perturbations. Then the LSR rate 
is easily related to the spectral density of the local mag- 
netic fluctuations, calculated at the NMR frequency de- 
termined by the external fieldi^i^i Also curious is the 
way nuclear spin diffusion proceeds in our system. The 
well-known treatment of nuclear relaxation by coupling 
to paramagnetic impurities plus nuclear spin diffusion — 
shows that there is a "spin diffusion barrier radius" be- 
low which neighboring nuclear spins cannot exchange en- 
ergy because the large dipolar field from the impurity 
brings them out of resonance. Here, instead, there is no 
such minimum radius for spin diffusion because nuclei at 
equivalent sites of different molecules are also magneti- 
cally equivalent (provided both molecules have the same 
electron spin orientation). 
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A. Spin-bath analysis and tunneling rate 

To apply the spin bath theory to the ^^Mn NMR in 
Mni2-ac, we begin by truncating the giant spin Hamilto- 
nian of the cluster to its tunneling-split ground doublet, 
and by taking as a basis for its subspace the m — ztS* 
projections of S along the z-axis, denoted by | fl-), | -IJ-). 
This restriction will be relaxed to consider higher ex- 
cited electron spin doublets when discussing thermally- 
assisted tunneling. Further, we assume that each central 
spin is coupled to N nuclear spins {Ife}, k = 1...N. 
The strength of each coupling is given by the quantities 
and fiiu^, which represent the part of the hyperfine 
coupling that does or does not change upon flipping the 
central spin, respectively (Fig. [T5)) . For nuclei in Mn*^^^ 
sites of Mni2-ac the hyperfine field Bhyp is exactly paral- 
lel or antiparallel to the direction of the cluster's z axis. 



so uj-^ = and cj 



= 7A'Bhyp. In Mn(2) and Mn^^) sites 
there's a small nonzero value of cj^ due to the orbital 
contribution to the hyperfine fieldi^ Conversely, for nu- 
clei such as ^H, which are subject to the vector sum of 
the dipolar fields from several surrounding clusters, we 
may expect and wj! to have comparable values. Let 
us define for each nuclear spin a number mt representing 
the spin projection of 1^ along the direction of the local 
hyperfine field Bhyp,fe. For nuclei irik = ±1/2, while 
for ^^Mn TOfc = —5/2 . . . + 5/2. Then the total hyperfine 

bias on the cluster is ^jv ~ ~'^^J2k=i '^fcwf . With this 
definition, < when the majority of nuclear spins 
is parallel to the local Bhyp,fc, thereby lowering the total 
energy of the system. Notice that, for a given orientation 
of the nuclear spins, changes sign whenever the elec- 
tron spin flips, since the direction of Bhyp does. Thus we 
define an absolute index of nuclear polarization in each 
cluster as 7-" = Cs^f^mk, with Cs = +1 when S is in 
the I ft) state, and C5 = — 1 otherwise. Each possible 
value of V deflnes a "polarization group" , and is inde- 
pendent of the electron spin state. Since the individual 
hyperflne couplings vary over a broad range (from ~ 1 
MHz for distant protons to 365 MHz in Mn^^^), the pos- 
sible values of the bias for each V are also widely 
spread, yielding a set of largely overlapping polarization 
groups. Globally, we may describe the coupled "central 
spin -I- spin bath" system by two manifolds of states, one 
for each electron spin state | ff) , | i^), split by hyperflne 
interactions into a dense band of states indexed by the 
nuclear polarization V, as shown in Fig. [TBI Calling 
'Pmax the maximum value assumed by V, Vmax — N if 
Ik — 1/2 Vfc. The proflle of the hyperflne bias distribu- 
tion can be calculated with the knowledge of the individ- 
ual couplings, and is well described by a Gaussian with 
half-width Eo = Ekiilk + l)/3/fc](cc;f/fc)2 ~ 0.082 

In addition to the hyperflne couplings, the S spins are 
also mutually coupled by dipolar interactions, which yield 
an additional bias — 2(7 /ibS • Bdip. The dipolar bias 
can be considered quasi-static in the sense that it remains 
essentially constant over time intervals that are long com- 




FIG. 15: (Color online) Scheme of the relative orientations of 



the hyperfine fields before (B^^^p j. 



) and after (B^^^p j.) the elec- 



tron spin flip, and the components of the hyperfine coupling 
that change (tjj!) or stay unchanged (u)^) at each tunneling 
event. The angle f3k is involved in the definition of k, the 
number of nuclei coflipping by "orthogonality blocking" , Eq. 
(Ei). 



pared to the typical timescale for the hyperflne bias fluc- 
tuations. The distribution of dipolar biases depends on 
the total magnetization of the sample and, in general, on 
its shape. For a demagnetized, ZFC sample of Mni2-ac, 
the dipolar bias distribution is described by a Gaussian 
with half-width E-d ~ 0.32 Finally, one may in gen- 
eral apply a static external fleld, Bz, along the z-axis, 
which produces an additional bias = ig^sSzBz- For 
zero external field and some typical nonzero value of 
the energy level scheme of a Mni2-ac cluster coupled to 
its nuclear spins would resemble the sketch shown in Fig. 

To analyze the behavior of this system with respect 
to incoherent tunneling of the electron spin, the crucial 
question to be answered is what happens to the nuclear 
spins when S suddenly changes direction. How many of 
the {Ifc} coflip with S? As extensively discussed in the 
PS literature, ^'^'^ there are two mechanisms by which 
nuclear spins may be fiipped by a tunneling event. First, 
a nuclear spin may coflip with S if the local hyperflne 
fleld does not exactly reverse its direction after S has 
tunnelled, since it would then start to precess around a 
different axis, hence the name "orthogonality blocking" 
or "precessional decoherence" for this mechanism. The 
number of spins coflipped this way is k, deflned as (see 
Fig.[l5]for/3fe): 



= Y\ cos Pk 



cos(2/3fc) 



B 



(1) 



B 



(2) 



hyp,fc hyp,fc 



IB 



(1) 



B 



(2) 



(9a) 
(9b) 



hyp,fc 1 1 hyp,/c 



The cosine factors in Eq. (|9ap . which are multiplied over 
all the bath spins, are the overlap matrix elements be- 
tween the initial and flnal bath states, i.e. {i\Uk\f), where 
Uk is the rotation operator of the k-ih. bath spin (Ref. 0, 
Appendix A. 2). Clearly, n depends only on the direc- 
tion of the hyperflne flelds, and not on the timescale of 
the electron spin flip. The nuclei in Mn^^^ sites do not 
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"low" dipolar energy 

FIG. 16: (Color online) Sketch of the hyperfine-split manifolds representing the energy of the m — ±S electron spin levels 
coupled to the nuclear spin bath. 



contribute to k since the Bhyp before and after the flip 
are exactly antiparallel, i.e., lo -L(Mn(^)) = 0. Conversely, 
nuclei in the ligands may give a large contribution 
because they are subject to the vector sum of the dipo- 
lar fields from several molecules, which does not entirely 
reverse direction when just one molecule flips. 

The other possibility is that the nuclear spins follow 
adiabatically the rotation of S. For this to happen, the 
"bounce frequency" of S, flo, has to be small or compa- 
rable with the nuclear Larmor frequencies. f2o is given 
here by the energy difference between the |m| = S' and 
1 771 1 = S' — 1 cluster spin states: since we are interested 
in FRMs, knowing that the resonance between m = —S 
and 771 = 8-1 states occurs at ~ 0.39 T (Ref. IH) 
yields HHq ~ 10 K, i.e., several orders of magnitude larger 
than {ujk}. Therefore, the nuclei cannot adiabatically fol- 
low the dynamics of S and the number of spins coflipped 
by this mechanism. A, is essentially zero. As a matter 
of nomenclature, this mechanism leads to what is called 
"topological decoherence" because the topological phase 
of the {Ik} becomes entangled with that of S.^'^'^ 

Combining the two flipping mechanisms deflnes a pa- 
rameter ^0 oc A-f K, which expresses how much the nuclear 
polarizations before and after the electron spin flip may 
differ for the flip to be likely to occur. Two opposite 
situations are sketched in Fig. [171 where we call V'^^^ 
and P'^' the nuclear polarizations before and after the 
electron spin flip, respectively. In any case the system 
has to tunnel between states at the exact resonance, but 
in case (a) the electron flip does not require any nuclear 
coflip (7-''^'' = 'P^'^^), while case (b) requires all nuclei 
to coflip {V'^^^ = — T'^^^), which is extremely unlikely. 
As a result, the expression for the tunneling rate con- 
tains a factor exp(— ^/^o) that describes precisely this 
restriction. From the above discussion it is clear that - 
at least in the absence of external transverse fields^^ - the 
main contribution to comes from nuclei, and that 




FIG. 17: Sketch of two resonant tunneling processes, differing 
in the number of required nuclear coflips. (a) The nuclear 
polarization is the same before and after the electron spin flip: 
this process has maximum likelihood, (b) All nuclei need to 
reverse their spin to conserve the total energy: this process is 
extremely unlikely. 



Co<{So,Sd}. 

In the presence of a dipolar bias, the tunneling tran- 
sition with highest probability, i.e. no coflipping nuclei, 
occurs when = ^^r (Fig. (TO)) . This means that a 
tunneling event effectively entails an exchange of dipolar 
and hyperfine energy. We may then distinguish between 
transitions that increase the hyperfine energy ("left to 
right" in Fig. I16p . occurring at a rate W|, and transi- 
tions that decrease it ("right to left " in Fig. fTB]) at a 
rate wi. The total tunneling rate is F''^ = {w[ + wi)/2. 
The PS expressions for F-^, generalized to the m-th elec- 
tronic doublet, given a (dipolar) longitudinal bias ^ on 
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the ground doublet, are)^ 
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A2 , 



where is the tunnehng matrix element of the m-th 
electron spin doublet, to be calculated by exact diagonal- 
ization of the giant spin Hamiltonian. The factor cj^' 
expresses the fact that the spin-bath mediated tunneling 
rate vanishes when A.^ 3> Eq"^\ i.e. when the spread of 
nuclear energies is not sufficient to sweep the hyperfine 
bias through the tunneling resonance. The parameters 
Eq"^\^^™^ are generalizations to arbitrary electron spin 
doublets of the quantities Eq and ^ defined before for 
the ground doublet |m| = S, while ^0 is assumed m- 
independent. To obtain the total tunneling rate through 
the m-th doublet we average Eq. pU|) over the distribu- 
tion of dipolar biases: 



1 



exp(-eV2i?lr^ 



2ttE 



(m) 



(11) 



In the real situation considered here, the spread of dipolar 
biases in the sample is much larger than the tunneling 
window allowed by hyperfine couplings, Ed ^ ^o- This 
means we can estimate by calculating the fraction 
of molecules with bias — f'™-' ^ c ^ 
may approximate r^(^) 
contribution of the molecules whose bias is larger than 
^0 and which tunnel at an exponentially small rate: 



Co"^ < e < for which we 

~ r^(0), and by neglecting the 
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Finally, the global spin-bath driven tunneling rate, , 
is obtained by summing over the m electronic doublets 
weighed with the appropriate Boltzmann occupation fac- 
tor: 



r^(T)c.l^ 



cxp 



P \ 9A2 



(m) ' 



(14) 



~£.o < £. < 'Co J we have c>c oc ^o- This immediately 
explains why the isotopic substitution of for yields 
a decrease in tunneling rate [Fig. [DJa)], since these are 
the nuclei that mostly contribute to ^o- 



B. Phonon-induced tunneling rate 

For comparison, we also discuss the case where the 
electron spin tunneling is caused by spin-phonon cou- 
plings. The phonon-driven tunneling rate though the m- 
th doublet, F^^, is related to the (T-dependent) broaden- 



ing of the electron spin states, Wm{T), by 



.70.71 



rt(0 



(15) 



The phonon-induced broadenings are obtained as a func- 
tion of the sample density, p, the sound velocity, Cg, and 
the uniaxial anisotropy parameter, D, asi^ 
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with s±i = (5 =F m){S ± m + 1)(2to ± 1)^ and s±2 — 
{S T m){S ±m + 1){S Tm- 1){S ±m + 2). 

Again, we calculate the fraction of molecules with high- 
est tunneling rate, a;^, as those whose bias is within the 
width of the Lorentzian function (jlSp : 



(17) 



and weigh the contribution of the m-th levels with their 
Boltzmann factor to obtain the total phonon-driven tun- 
neling rate: 



F^r)c.i^exp 



k^T) ^^A^^ + h^wUT)-'^ ' 



Contrary to the nuclear-driven case, here each individual 
m-th doublet tunneling rate F^ is T-dependent by itself, 
besides being weighed by Boltzmann factors. This means 
that the phonon-driven tunneling rate never shows a T- 
independent plateau, even at very low-T when tunneling 
occurs only through the ground doublet. 



where {Em} arc the average energies of the m-th doublets 
and Z is the partition function. Notice that the spin-bath 
driven tunneling rates are individually T-independent: 
the temperature enters only in the Boltzmann factors 
for the occupation of the m-th doublets, and thereby in 
their contribution to the global tunneling rate F^(r). 
Also, since Pm{£,) is essentially constant in the interval 



C. Tunneling rate 



longitudinal spin relaxation 
rate 



To relate the electron spin tunneling rates, F^^'^(T), to 
the observed LSR rate, we apply the remarks made 
above on the behavior of the nuclear spins upon a sudden 
change of the electron spin direction. In particular, we 
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shall compare theory and experiments for the Mn^^^ site, 
where (ideally) the hyperfine coupling would be strictly 
scalar.— This implies that the ^^Mn spin at Mn^^^ sites do 
not coflip by "precessional decoherence" , since uj-j: ~ 0, 
neither do they coflip by "topological decoherence" due 
to the very small values of a;|/rio- Thus, each electron 
spin tunneling event corresponds to the inversion of the 
populations of the local nuclear Zeeman levels. The LSR 
rate arising from this situation is easily obtained and can 
be found in the literature^^i^i^i^^ but we repeat here the 
derivation because it will allow us to point out exactly 
why no known theory can explain our data. The answer 
is most easily obtained for nuclear spins / = 1/2, but 
remains valid for arbitrary spin values. 

We start by writing a master equation for the popu- 
lations of the nuclear Zeeman levels relative to the local 
hyperfine field direction^ calling N-^- the number of nu- 
clei in the excited Zeeman state and 7V_ those in the 
ground state. For simplicity, since the internal equilib- 
rium is reestablished within a time T2 ^ tt after each 
tunneling event, we assume that just before tunneling all 
clusters have the same values of and iV„ , neglecting 
fluctuations around the mean values. If a fast-relaxing 
molecule tunnels at time t, the polarization of its own 
nuclei is abruptly inverted. Each time a tunneling tran- 
sition lowers the energy of the local nuclei, which occurs 
at a rate wi, then N- nuclei have been added to the to- 
tal number of nuclei in the Zeeman ground state. After 
a time T2 this decrease in local hyperfine bias has been 
redistributed over the sample: calling a^pRM the fraction 
of FRMs over the total, then the tunneling event has in- 
creased to N+ + x-pkmN^. The same reasoning holds 
for transitions that increase the hyperfine bias. The mas- 
ter equation is therefore: 



diy+ 
dt 



(19) 



From here the LSR rate can be obtained by standard 
textbook calculationsi^i Writing = {N + n)/2 and 
N_ = {N - n)/2, Eq. ^ becomes: 



= XYYiMN{wi - W|) - XYKMn{wi +w^), (20) 



which can be rewritten as: 



An 

— = 2W{no - n) 



W = XFRM 



N,<j> 



(21) 
(22) 

(23) 



where no is the equilibrium nuclear polarization and W is 
the desired LSR rate, since the solution of (PT|) is precisely 
of the form n{t) = n(0) - [n(0) - no][l - exp(-2W^t)]. 

We now attempt to fit the measured LSR rate in zero 
field at the Mn'^^ site in ZFC sample, using Eqs. dH]), 



(IT5)) and ((23)) . To this end we calculate the energy levels 
scheme of the FRMs using the effective spin Hamiltonian: 

Hfrm = -DSl + E{Sl - SD - C{SX + St). (24) 

Unfortunately, very little is known about the parameter 
values in ([24|l . To the best of our knowledge, it's not even 
established whether FRMs in Mni2-ac have lowest total 
spin state S* = 10 or, for instance, S* = 9. The analysis of 
a Mni2 variant containing only FRMs^^ seemed to sup- 
port an S" = 10 ground state, but it's not clear to what 
extent the FRMs in Mni2-ac have the same properties as 
those analyzed in Ref. [t^- For instance, Ref. 77 finds 
the first level crossing transition, i.e. the value of longi- 
tudinal field at which Es = E_s+i, at B^-^+'^ ~ 0.27 
T, quite different from the value observed for the actual 
FRMs in Mni2-ac, Bf'-^+i ~ 0.39 T.^s We shaU try 
both 5 = 10 and S = 9 and discuss how the different 
behaviors compare to the experimental data. To avoid 
having too many fitting parameters, we choose to keep C 
fixed at the value commonly used for the majority species 
of Mni2-ac, C = 4.4 x 10~^ Ki^ The uniaxial anisotropy, 
D, is obtained by imposing the condition B^ ~^'^^ ~ 0.39 
T,'^^ yielding D = gfXB x 0.39 = 0.524 K. Ahhough we 
tried adding also a fourth order term, —BS*, it turned 
out that the best fits are obtained by leaving B = 0, so 
we shall not discuss this further. The rhombic anisotropy 
term, E, is used as the actual fitting parameter since it 
most directly influences the value of the tunneling split- 
tings 2 Am and thereby the tunneling rates F^'*^. Car- 
retta et alJ^ showed that the effective A is extremely 
sensitive to the gap between the lowest lying total spin 
manifolds ("S-mixing"), and could explain why the ob- 
served Landau-Zener tunneling probabilities in Feg are 
much larger than what would be expected on basis of the 
spin Hamiltonian parameters for the S" = 10 manifold.— 
A small gap between lowest lying total spin manifolds is 
quite expectable for FRMs, so the values of E we need to 
justify the F^'"^ extracted from experiment should not be 
taken literally as an estimate of the anisotropy parame- 
ter. In other words, the values of E used in our calcu- 
lations account also for the possible S-mixing due to an 
energetically close manifold with different total spin S, 
and do not necessarily correspond to the values that one 
would obtain from neutrons scattering or EPR experi- 
ments. We take as fixed parameters the sound velocity 
Cs = 1.5 X 10^ m/s (Ref. 80), the density p = 1.83 x 10^ 
g/m3 (Ref. (HI), Ed = 0.32 K (Ref. M) and Eq = 0.082 
K (Ref. ItoI ). whereas ^0 is allowed to vary. When com- 
paring nuclear- and phonon- driven tunneling rates, we 
impose the same parameters for the spin Hamiltonian 
([24| . The results of the calculations are shown in Fig. 
[T51 for the set of parameters given in Table |TT1 

By looking at the theoretical curves alone we find that, 
given a fixed set of parameters for the FRMs' spin Hamil- 
tonian (|24[) . the nuclear-driven tunneling process always 
dominates over the phonon-driven one, both in the low-T 
and in the high-T regime. The situation may be reversed 
in the high-T regime by assuming the sound velocity is 
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FIG. 18: (Color online) Calculated nuclear spin-lattice relaxation rates, W , as a function of temperature, for spin-bath [panels 
(a) and (c)] and phonon [panels (b) and (d)] mediated tunneling. In panels (a) and (b) a total spin 5 = 10 is assumed for the 
FRMs; in panels (c) and (d), 5* = 9. The complete parameter sets are given in table |IT1 Black dots: experimental data; thick 
black lines: calculated W\ thin lines: contributions of the m-th electron spin doublets to the total rate. 



TABLE II: Parameter values for the calculations shown in 
Fig. 1181 In bold are given the free fitting parameters. The 
tunneling matrix element in the lowest doublet, As, is ob- 
tained from the Hamiltonian (|24p. i.e. is not introduced by 
hand. 
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0.204 
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0.524 
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0.524 
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10 
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0.524 
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lower than the literature value used here, but Y"^ > 
would never hold at low-T under realistic circumstances. 
A comparison with the experimental data shows that 



both nuclear- and phonon- driven mechanisms yield a 
correct slope of W{T) in the thermally- assisted regime, 
T > 0.8 K, whereas the phonon process can never re- 
produce the low-T plateau. An almost perfect fit of the 
data is obtained by assuming that the FRMs have a to- 
tal spin S = 9, while the S' = 10 case has a T-dependent 
region systematically starting at too high temperatures. 
In the nuclear-driven case, we used an optimal value of 
^0 — 10 mK which seems very reasonable, since the main 
contribution to arises from the coupling to protons. 
Conversely, the values used for E appear very high, since 
to fit the NMR data we need to assume D/E ~ 2.5. 
As mentioned before, however, such a high value of E 
should not be interpreted as the spectroscopic rhombic 
term in the spin Hamiltonian, since it is used here as 
the parameter that tunes the tunneling splitting of the 
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ground doublet, and therefore incorporates the effect of 
S'-mixinglS if a manifold with different total spin value is 
energetically close to the S ground state. 

An important remark is that our calculations, which 
only account for tunneling fluctuations as a source of nu- 
clear LSR, can accurately reproduce the observed LSR 
rate also in the thermally activated regime, T > 0.8 K. All 
the previous NMR experiments in that regime^iiSiii^iSi 
have been interpreted in terms of the "intrawell" electron 
spin fluctuations, that arise from thermal excitation of 
the electron spin state on the same side of the anisotropy 
barrier [see Fig. mb)]. As pointed out by Goto et al.^ 
that reasoning is inappropriate when applied to the LSR 
of nuclei belonging to Mn'^+ ions, since the hyperfine cou- 
pling tensor is diagonal. In that fluctuation in the 
z projection of the electron spin does not result in a fluc- 
tuating field perpendicular to the nuclear quantization 
axis (which is z itself), and cannot account for longitu- 
dinal nuclear spin relaxation. As we show in Fig. 1181 
including the effect of electron spin tunneling through 
|m| < S doublets solves what appeared to be a para- 
dox, since the tunneling fluctuations induce nuclear LSR 
also in the absence of non-diagonal hyperfine coupling 
terms. For the nuclei in Mn'^^ ions, the intrawell elec- 
tron spin transitions do provide an additional channel for 
nuclear relaxation, which may explain why the LSR rate 
at T > 1 K in the Mn'^+ sites is larger than that in the 
Mn4+ onesi^^ 



D. Thermal equilibrium 

In the preceding discussion, it may seem that we have 
not explicitly used the condition that the nuclear spins 
are in thermal equilibrium, which is what we observe in 
the experiment. This condition, however, is automati- 
cally implied in the application of Eq. (I23p to the LSR 
rate. For Eq. ([23)) to actually represent the rate at which 
the nuclear spins exchange energy with a thermal bath 
and thereby return to the equilibrium magnetization af- 
ter a perturbing NMR pulse, one needs to include the 
detailed balance condition: 

— =exp -— — . (25) 

In other words, if the nuclear spin temperature has to 
reach equilibrium with the thermal (phonon) bath via 
the process described by the rates and the latter 
must satisfy (f25|) . The detailed balance condition is of- 
ten taken for granted, but in this case one needs to be 
more careful, and wi represent the rates of electron 
spin transitions that increase or lower the nuclear energy, 
respectively. The crucial point is that both are rates for 
tunneling transitions, which occur when the total energy 
of the "electron plus nuclear spins" system is the same 
before and after the electron spin flip. Thus, the differ- 
ence between w-f and is simply that, e.g., is the 
rate for a tunneling transition that increases the nuclear 



spin energy while reducing the electronic one. That is, 
after the flip most of the nuclei are oriented against their 
local hyperfine field, while the electron spin is favorably 
aligned with respect to the local field (in particular the 
dipolar one, when Bz =0). W| does the opposite and, in- 
terestingly, this means that the instantaneous local spin 
temperature (local referring to the nuclei belonging to a 
specific molecule that has just flipped) is negative. This 
situation is clearly very different from the standard NMR 
picture of nuclear relaxation by coupling to paramag- 
netic centers, where the latter make spin-phonon tran- 
sitions between Zeeman-split levels having different ther- 
mal populations. 

Now we can summarize the meaning of our experimen- 
tal results for the description of electron spin tunneling 
in the presence of a nuclear spin bath: 

(i) The Prokof'ev-Stamp theory of the spin bath, as 
developed so far and reviewed in i ^V Al quantitatively 
and qualitatively reproduces the nuclear LSR rate in the 
whole temperature regime of our measurements, by as- 
suming that the LSR is due to tunneling events in a mi- 
nority of fast-tunneling molecules; 

(ii) The additional observation that the nuclear spins 
are in thermal equilibrium with the phonon lattice at all 
temperatures implies that the rates and wi must he 
different. For this to happens, it is necessary to explic- 
itly include the role of spin-phonon interactions in the 
nuclear-spin mediated tunneling process. Importantly, 
the results of the calculations shown in Fig. [TBTb.d) in- 
dicate that it is not sufficient to attribute the thermal 
relaxation to a phonon-assisted tunneling process as de- 
scribed in WBl working "in parallel" to the nuclear-spin 
mediated tunneling process. At the lowest temperatures, 
even the longest thermalization times observed in our ex- 
periments (i ^IV Ap arc still much shorter than (F*^)"^ as 
calculated from phonon-assisted tunneling alone ( WBp . 
thus reinforcing the need for a theory that includes 
nuclear-spin and phonon mediated tunneling at the same 
time] 

(iii) Our statement that the nuclear relaxation has 
to be mediated by inelastic electron spin tunneling pro- 
cesses, is further supported by specific heat experiments 
that show that a system of dipolarly-coupled tunneling 
molecules can relax to the long-range ferromagnetically 
ordered state, provided that the tunneling rate is fast 
enough for the experimental detection of the ordering 
anomaly. This means that there is a mechanism for 
the ensemble of electron spins to find its thermodynamic 
ground state, even at temperatures so low that the relax- 
ation can only proceed by quantum tunneling. We argue 
here that such inelastic tunneling mechanism is the same 
mechanism that is responsible for the thermalization of 
the nuclear spins. Indeed, by extending the specific heat 
measurements below the ordering temperature, one even 
observes the equilibrium specific heat contribution of the 
nuclear spins, meaning that both the electron spins and 
the nuclear spins can attain thermal equilibrium within 
the time-scale of the specific heat experiment (10 - 100 
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s) at all temperatures reached. 

All the work presented here is dedicated to the small- 
Aq, incoherent tunneling regime for the central spin. 
Having shown that the description of the nuclear-spin 
mediated tunneling is incomplete without the inclusion 
of spin-phonon couplings, some concerns may be raised 
also on the current description of the spin bath effects 
on the electron spin in the large- Aq, coherent tunneling 
regime This work cannot address that issue, and 
we think that the answer will have to come from low-T 
NMR experiments in large transverse field, and quanti- 
tative analysis of pulscd-ESR experiments. 

VI. CONCLUDING REMARKS 

The purpose of the research presented here is to il- 
lustrate and analyze a prototypical example of quantum 
tunneling of a macroscopic variable (the giant spin of 
single-molecule magnet) in the presence of a spin bath 
environment. Instead of looking at the macroscopic vari- 
able itself and deducing the effect of the environment on 
its dynamics, as is most often done, we have directly ob- 
served the behavior of the spin bath by means of low-T 
NMR experiments. 

We have provided compelling evidence that the lon- 
gitudinal nuclear spin relaxation in the ^^Mn nuclei of 
Mni2-ac is driven by electron-spin quantum tunneling 
fluctuations. The nuclear LSR rate, W, indeed con- 
tains all the features that are expected to be associated 
with tunneling of the molecular spin: i) A T-independent 
plateau of the LSR rate for T < 0.8 K; ii) A strong depen- 
dence of on a longitudinal magnetic field, that destroys 
the resonance condition for electron spin tunneling; iii) 
The slowing down of the nuclear LSR upon isotropic sub- 
stitution of by in the ligands, by an amount iden- 
tical to the slowing down of the quantum relaxation of 
the magnetization observed in similar systems. Because 
of the short timescale of the observed LSR, we argued 
that the tunneling fluctuations must take place in a mi- 
nority of fast-relaxing molecules, which are indeed known 
to be present in Mni2-ac. For these fluctuations to re- 
lax the nuclear magnetization in the entire sample, an 
additional mechanism is required which equilibrates the 
nuclear spin polarization across neighboring molecules, 
i.e. intercluster nuclear spin diffusion. Our data on the 
transverse nuclear spin relaxation show that the inter- 
cluster spin diffusion is indeed present and effective. All 



the above observations confirm and support the picture 
of nuclear-driven quantum tunneling of magnetization as 
originally formulated by Prokof 'ev and Stamp. However, 
a crucial outcome of our experiments is the demonstra- 
tion that the nuclear spins are in thermal equilibrium 
with the lattice phonons down to the lowest tempera- 
tures, where only quantum tunneling fluctuations of the 
electron spins are still present. This observation cannot 
be explained within the present theory of the spin bath. 

The implications of our results are potentially very pro- 
found, particularly because of the growing interest to- 
ward a coherent manipulation of spins for quantum in- 
formation processing. The spin-bath environment, de- 
scribing localized two-level systems, has been repeat- 
edly identified as the most important source of decoher- 
ence in solid-state qubits. This includes superconducting 
systems,— quantum dots,— NV centers in diamondSl 
and, of course, molecular magnetsi^ We have investi- 
gated here the incoherent tunneling regime, but the the- 
oretical formalism to describe the coupling between cen- 
tral spin and spin bath is identical in the case of coherent 
spin dynamics. Therefore, the main finding of our work 
- that the role of phonons in the nuclear-spin mediated 
tunneling is currently lacking a proper description - sug- 
gests that also the contribution of the nuclear spin bath 
to the decoherence rate of realistic spin qubits may need 
to be revisited. 
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